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ABSTRACT 

Context. Classical Be stars are rapid rotators of spectral type late O to early A and luminosity 
class V-III, wich exhibit Balmer emission lines and often a near infrared excess originating in 
an equatorially concentrated circumstellar envelope, both produced by sporadic mass ejection 
episodes. The causes of the abnormal mass loss (the so-called Be phenomenon) are as yet un- 
known. 

Aims. For the first time, we can now study in detail Be stars outside the Earth's atmosphere with 
sufficient temporal resolution. We investigate the variability of the Be Star corot-id 102761769 
observed with the CoRoT satellite in the exoplanet field during the initial run. 
Methods. One low-resolution spectrum of the star was obtained with the INT telescope at the 
Observatorio del Roque de los Muchachos. A time series analysis was performed using both 
cleanest and singular spectrum analysis algorithms to the CoRoT light curve. To identify the 
pulsation modes of the observed frequencies, we computed a set of models representative of 
corot-id 1027617 39 by varying its main physical parameters inside the uncertainties discussed. 
Results. We found two close frequencies related to the star. They are 2.465 c d~' (28.5 yuHz) and 
2.441 cd~' (28.2 yuHz). The precision to which those frequencies were found is 0.018 cd~' (0.2 
yuHz). The projected stellar rotation was estimated to be 120 kms~' from the Fourier transform 
of spectral lines. If corot-id 10276 17 39 is a typical Galactic Be star it rotates near the critical 
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velocity. The critical rotation frequency of a typical B5-6 star is about 3.5 c d~'(40.5 /jHz), which 
implies that the above frequencies are really caused by stellar pulsations rather than star's rota- 
tion. 

Key words. Stars: early-type - Stars: emission-line, Be - Stars: individual: corot-id 10276 17 39 
- Stars: rotation - Stars: oscillations 



1. Introduction 



Classical Be stars are rapid rotators of spectral type late O to early A and luminosity class V-III, 
wich exhibit Balmer emission lines and often a near infrared excess originating in an equatorially 
concentrated circumstellar envelope, both produced by sporadic mass ejection episodes. The causes 
of the a bnormal mass loss (t he so-called Be phenomenon) are as yet unknown. In spite of their high 



V sin i dFremat et al. 



2005I) . rapid rotation alone ca nnot explain the ejection episodes as most Be 

for a comprehensive 



Porter & Rivinius 2003 



stars do not rotate at their critical rotation rates (see 
review of Be stars). 

High-resolution, high signal-to-noise spectroscopic observations have been analyzed to demon 



strate short-term variations are rather common among early -type Be stars (IBaade 



Porri & Stalio 



1988 



Roquet et al. 



1992; 



Neiner et al. 



1982, 



1984 



2002). The observed line profile variabil- 



ity (LPV) is characterized by moving bumps traveling from blue to red across the line profile on 
timescales ranging from min utes to a few days. T he phenomenon has also been observed in O 
stars and 6 Set variables (e.g. 



Kennelv et al. 



19921) . among others. Non-rad ial pulsations (NRP) 



have been proposed as an ex planation of the LPV observed in hot stars (e.g. 



Smith 



1977 



Baade 



1982; 



Vogt & Penrod II 19831) . NRP produce LPV thanks to the combination of the Doppler dis- 



placement of stellar surface elements with their associated temperature variations due to the com- 
pression/expansion caused by the passage of waves through the photosphere. NRP could be the 
additional mechanism required for a rapidly rotating B st ar to become a Be s tar, that is to trigger 
the Be phenomenon by means of mass ejection. Indeed. iRivinius et al. I dl998) found multiperiod- 



icity in yu Cen and showed that a correlation exists between mass ejection episodes and the beating 
pattern of the multiperiodicity. The finding of new cases of Be stars for which beating periods of 
multiperiodic NRPs coincide with matter ejections would help us to confirm this model. The peri - 



odic variability of the star has been reproduced in detail by NRP modeling (IRivinius et al. 



2001). 



Short-periodic LPV of other Be stars have also been modeled using NRP. Observations with the 
MOST satellit e showed that multiperiodicity due to NRP is a rather common phenomenon among 
Be stars (e.g. IWalker et al. 1 12005c iDziembowski et al. 112007b . Nevertheless, to measure frequen- 



cies with great accuracy in stars requires both high photometric precision and high time-frequency 
resolution. The CoRoT (Convection, Rotation and planetary Transits) satellite opens for th e first 



time in history the possibility of fulfilling this goal. CoRoT dBoisnard & Auvergne 



2004) is an 



experiment of astronomy dedicated to seismology and the detection of extrasolar planet transits. 
It was launched on December 2006 in an inertial polar orbit at an altitude of 897 km. The instru- 
ment is fed by a <p — 27cm telescope. Its scientific program is three-fold consisting of: (1) The 
seismology core program (SISMO), which concerns the seismic study of -10 bright (6 < V < 9) 



* The CoRoT space mission was developed and is operated by the French space agency CNES, with partic- 
ipation of ESA's RSSD and Science Programmes, Austria, Belgium, Brazil, Germany, and Spain 
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well-selected pulsating stars per observed field. Those measurements are sampled with a cadence 
of 32 seconds, which allows very accurate frequency analysis. (2) The exoplanet core program 
(EXO), which concerns the search for exoplanets around ~ 12000 faint (11.5 < V < 16) stars per 
field. The sampling of the EXO data has a cadence of 512 seconds. (3) The additional program 
(AP), which concerns the study of a few hundreds faint (11.5 < V < 16) variable stars per EXO 
field. Within the core exoplanet data, there are hot stars among which there will certainly be new 
pulsators, enriching the statistics of the NRP phenomenon. In this work we report on the analysis of 
the star corot-id 102761769. Because of a certain instability of the attitude control system (satellite 
jitter), the analysis of some of CoRoT stars requires a careful treatment. In the next section, we 
present the light curve as coming from the reduction pipeline and explain the measured instrumen- 
tal noise. In Sect. 3, we discuss the time analysis of the light curve and how the information was 
extracted. In Sect. 4, we describe the spectroscopy ground-based observations of the star and an 
estimate of the rotation period. In Sect. 5, we discuss the frequencies found and compare them with 
the critical rotation velocity. In Sect. 6, the conclusions of this work are presented. 



2. CoRoT photometric observations 



The data were taken in the exo field during the initial run (IRaOl) during 54.6 days with an initial 
cadence of 5 12 seconds. An automatic procedure (alarm mode) triggers a higher cadence in the case 
of suspect detection of a planetary transit. After 3.13 days, CoRoT began to record the data of the ID 
102761769 star with a cadence of 32 seconds (Fig. 1-top). The coordinates (J2000.0) of the star are 
RA=06 /I 45'" 02M and DEC=-01° 50' 08". Q3 and V= 13.288 mag. The star is also cited in cat- 



alogu es mentioning as emission-line object ( Ko houtek & Wehmever 



1999; 



Robertson & Jordan 



1989). The light curve contains 129 305 photom etric observations. Looking for oscillations, we 
first applied the cleanest algorithm (IFoster Nl995l> to the data. Figure 2 shows the first cleanest 
power spectrum. We identified several frequencies, the most easily distinguishable ones being re- 
lated to the satellite's orbit and its harmo nics. The close distance from Earth causes perturbations 
on the satellite leading to satellite jitter dAuvergne et al. 1120091) . The main sources are eclipses, 
the variation in the Earth's gravity and magnetic field, South Atlantic anomaly, the Sun and Earth 
emissivity, Earth's albedo, and objects in low Earth orbit (LEO). Figure 3 shows a phase diagram 
of the raw data with the 6184 seconds CoRoT satellite orbit period. CoRoT pipeline have a flag 
for bad points. Some of the points were corrected by the pipeline, while others were eliminated, 
which is why we have missing points. The gaps in Fig. 3 were caused by points eliminated from 
the automatic pipeline reduction mainly by satellite jitter. The cluster of points in Fig. 3 are the re- 
maining points that the pipeline could not correct. We eliminated 480 points greater than 3 cr from 
the average (0.4% of the total points). We also eliminated the first 355 observations (0.3% of the 
total points) taken with a larger cadence and filled the gaps with the average value (Fig. 1 - middle). 



The lo wer frequencies shown in the power spectrum in Fig. 2 are due to hot pixels dAuvergne et al. 



2009). There is a prism mounted in front of the CoRoT exoplanet channel camera detectors, which 



allows CoRoT to gather 3-colour light curves. We analyzed the total sum of the light curves (the 
so-called white light curve) to highest quality relation in terms of signal/noise. Each color light 
curve falls on different pixels, what allows us to check for cosmic rays. The bump around date 
2600 (Fig. 1 - middle) is only found in the red channel, and is indicative of cosmic-ray source 
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noise. The lower frequencies were eliminated by a polynomial fit to the data and the missing points 
were interpolated (17715 in total). An average of each 100 points was finally evaluated by reusing 
each point 100 times and then selecting 1 point every 100 points for additional analysis (Fig. 1 - 
bottom). The total number of points in Fig. 1 (bottom) is 1466. The regulary spaced intervals of 
data series were essential to the time analysis described in the next section. 
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Fig. 1. Top:corot-id 102761769 star light curve. Middle: Data after outlier removal at three cr level. 
Bottom: Moving average every 100 points. 



3. Time analysis 



We applied the singular spectrum analysis (SSA) (IGhil et al. 1120021) to the data in Fig. 1 -bottom. 
SSA is designed to extract information from noisy and evenly spaced time series. The method pro- 
vides insight into the unknown or only partially known dynamics of the underlying system that 
generated the series. The procedure takes an univariate record and produces a multivariate set of 
observations by using lagged copies of a single time record, we are able to define the coordinates 
of the phase space that approximate the dynamics of the system. SSA has been applied extensively 
to the study of climate variability as well as other areas in physical and life sciences. To create 
the singular spectrum we have to sort the eigenvalues of the trajectory matrix. Each eigenvalue is 
associated with an oscillation in the time series. The full line in Fig. 4 (top) shows the singular 



spectrum applied to the Be star corot-id 102761769. This plot is the stren gth of the sorted eig enval 



ues. We performed a 100 runs of Monte Carlo SSA shuffling the data (lAllen & Smith 



1996). The 



dashed line in Fig. 4 represents the 3<x confidence limit for this analysis. This means that the first 
four eigenvalues, defined by our Monte Carlo analysis, almost entirely encompass the signal of our 
time series. We reconstructed the time series (filtered series) with this limit (3 rd row of Fig.|4|. For 
comparison, we plotted again the bottom curve of Fig.[TJon a similar scale as our filtered series in 
the 2 nd row of Fig. [4] The reconstructed signal is a typical portrait of a beating process as expected 
for a Be st ar. Since SSA provides only the most significant eigenvalues found into data, we applied 
cleanest dFoster 1 1 1 9951) to identify the equivalent frequencies. Cleanest projects the data into an 
orthogonal subspace of trial functions (sines, cosines, and a constant function), while the projec- 
tion provided by SSA is onto the data itself. Figure [4] (bottom) shows the cleanest spectrum of the 
reconstructed time series. The cleanest spectrum follow a chi-square distribution with r-1 degrees 
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of freedom, where r is the number of trial functions (IFoster 



1996). The frequencies at 2.465 c d 



and 2.441 cd~' are more than 99.9% above the significance limit. Figure [5] shows the phase dia- 
grams for the two frequencies found for both smoothed and filtered series. We had already found 
those frequencies in the original light curve. To ensure that these frequencies were not introduced 
by numerical manipulation, we searched for them in each step of our analysis. Figure [6] shows a 
zoom in the 2.5 c d -1 region in the cleanest spectrum. From top to bottom we display, the light 
curve (lc), the lc without orbit perturbations, the lc with interpolated gaps, and the reconstructed 
time series. For all the light curves, the signal is above the significance limit. The difference be- 
tween the two close frequencies is 0.024 cd l above the limit of the frequency resolution, which 
is 0.018 c d~'. The amplitude found for both frequencies is 0.14 + 0.06 mmag using the algorithm 
cleanest. The CoRoT photometric performan ce is 7.7xl0~ 4 magnitudes in lh as faint as V=15.5 



dAuvergne et al. 



2003). 



Aigrain et al. I (12009b evaluated the noise per exposure on 2h timescales 



for the IRaOl and found the photometric performance to be close to the pre-launch specification. 
They found an empirical relation between 2 h noise and R-magnitude. Applying this empirical 
relation for the corot-id 102761769 star and taking all data analyzed (53.47 days), we found that 
o~53Aidays =0.01 mmag. However this light curve is particulary noise because of satellite jitter. The 
standard deviation of the mean from our data is 0.01 mmag. The outliers in the same figure are still 
caused by the jitter that was not totally eliminated from our analysis. This is why it is difficult to 
see the 2 frequencies from the 2 nd row of Fig. [4] 
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Fig. 2. Cleanest spectrum of the light curve (Fig. 1 - top). Most of the peaks are harmonics of the 
satellite orbital period due to the gaps in phase with the orbit (see Fig. 3). Among the significant 
frequencies, only two of them are not related to the the satellite orbital period harmonics. 



4. Spectroscopy ground-based observation 

A spectroscopic observation of the corot-id 102761769 star was performed at Observatorio del 
Roque de los Muchachos on March 30 2008 with the Intermediate Dispersion Spectrograph (IDS) 
at the 2.5m Isaac Newton Telescope. The spectrum was taken with 300 seconds of exposure time 
using the 300 1/mm grating centred at 5460 A and covering the region from 3337A to 8366A 
and resolution of 2 A (Fig. [7J». Figure [8] shows views of the H-a, H-/3, and He 4471A/Mg 4481A 
regions. The reduction was accomplished using IRAF software. The Ha line is seen in emission 
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Fig. 3. Phase diagram of the satellite folded modulo the orbital period (102.9 min). The gaps that 
appear in phases 0.3 and 0.95 are data eliminated by the automatic pipeline reductions due mainly 
to satellite jitter. 
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Fig. 4. Top: The full line shows the singular spectrum and the dashed line shows an average of 100 
runs of Monte Carlo simulation. 2 nd row: Smoothed time series (same as Fig. 1 - bottom). 3 rd row: 
Reconstructed time series using only the leading 4 eigenvalues, showing a strong beating effect. 
Bottom: Fourier transform of the reconstructed series. Note that the frequencies around 2.4 cd~' 
are by far the dominant ones. 



and H/3 is partially filled with emission. The spectral characteristics are indicative of a spectral type 
B5-6 IV- Ve based on the EWs of Hel lines, H5, Hy, and Mg II (,14481 A). The CoRoT catalog 
indicates for the star an A0V spectral type and lu minosity class. This classification was obtained 
from photometric data using a statistical method (iDeleuil et al. 1120091) . Our spectral classification 
based on spectroscopic data provides a more precise classification. The usual caution must be taken 
concerning the spectral classification of Be stars, which can be influenced by the circumstellar 
envelope and the gravitational darkening. The rotation velocity was estimated by evaluating the 
first zero of the Fourier transform of several He lines. We conclude that the star rotates with V 
sin i = 120 km s . This is a maximum value for the projected stellar velocity, because the spectral 
resolution is about 2 A, which corresponds to * 120 km s _1 . 
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Fig. 5. Phase diagrams. From top to bottom: Smoothed data with frequency = 2.465 cd~', SSA 
result with frequency = 2.465 c d -1 , smoothed data with frequency = 2.441 c d~', SSA result with 
frequency = 2.441 cd -1 . 




Fig. 6. Cleanest Spectra between 2.4 and 2.5 cd from top to bottom: the light curve (lc), lc 
without orbit perturbations, lc with interpolated gaps, reconstructed time series. 
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Fig. 7. Spectroscopic ground-based observation of corot-id 102761769 star. 



5. Discussion 



The two dominant frequencies found in Sect. 3 are /i = 2.46 5 c d 1 and f% = 2.44 1 cd . Taking 



M« = 3.8M G and R* = 2.6 / ? for a non-flattened B 5.5V star dHuang & Gies 112006 ) and using the 
expression of IZorec et al. I (12005b . the critical rotation velocity is 528 kms 1 and v fn , = 4 cd -1 . 



8 Emilio et al.: Photometric variability of the Be star corot-id 102761769 




Fig. 8. Close-ups of interesting spectral regions. From top to bottom: Ha, H/3, He 447 1 A/Mg 448 1 
A 



Fremat et al. 



(2005) found that typical Galactic Be stars rotate on average at 88 % of their critical 
rate. If that is the case for corot-id 102761759, v rot - 3.5c d _I and frequencies near f\ and fa found 
in this paper are much lower than the probable rotation frequency, which implies that they are true 
stellar NRP pulsation frequencies. The determined V sin i and the above critical velocity set the star 
inclinatio n angle to be » 15 degrees. From the Ha Ew (6.6 A) and FWHM (6.2 A), the relations 
given by 



Dachs et al. 



( 1986) allow us to estimate the inner radius of the region emitting the bulk 
of the line to be R a « 4.4/?,. On the other hand, the absorption cusps seen at the central region of 
the H/3 line indicate that the inner radius of the H/3 line-forming region is Rr ~ 1 .95R*. 



6. Conclusions 

We have analyzed 54.6 days of photometric data of the corot-id 102761769 target (V = 13.08, 
B5-6 IV- Ve star) obtained during the IRaOl run. One spectroscopic observation was also per- 
formed at Roque de los Muchachos observatory (La Palma, Spain), covering the spectral region 
/U=3337A-8366A. The Ha line was observed in emission and H/3 was found to be partially filled 
with emission. The projected rotational velocity is V sin i — 120 km s _1 , as estimated by means of 
the Fourier transform of Hel line profiles. The time series analysis was based on both cleanest and 
SSA algorithms. We found two dominant frequencies, f\ = 2.465 cd~' and fa = 2.441 cd~'. A 
strong beating between these frequencies was uncovered when the data were reconstructed by de- 
tecting the four most significant eigenvalues found applying SSA. If the star rotates near the critical 
frequency, the frequencies f\ and found in this paper are much lower than the probable rotation 
frequency. 

Acknowledgements. This work is based on observations made with the INT Telescope under the Spanish Instituto de 
Astrofisica de Canarias CAT Service Time. We thank the support astronomers team at the IAC, and, in particular to Rafael 
Barrena, for making the observations and providing us all the information needed to work on the data. This work was 
supported by FAPESP (Grant 05/58746-1 and 06/50341-5) and CNPq (Grant 477.813/2007-0 and 480.658/2008-0). 



References 

Aigrain, S., Pont, F., Fressin, F, et al., 2009, A&A, 506, 425 
Allen, M.R., & Smith, LA. 1996, Journal of Climate, 9, 3373 

Auvergne, M., Boisnard, L., Buey, J.-T.M., et al., 2003, in Proc. SPIE 4854, ed. J. C. Blades, & O. H. W. Siegmund, 170 
Auvergne, M., Boisnard, L., Lam-Trong, T, et al., 2009, A&A, 506, 411 
Baade, D. 1982, A&A, 105, 65 
Baade, D. 1984, A&A, 135, 101 

Boisnard, L., & Auvergne, M. 2004, CoRoT mission engineering, IAC-04-IAFQ. 1.01 



Emilio et al.: Photometric variability of the Be star corot-id 102761769 



9 



Dachs, J., Hanuschik, R., Kaiser, D., Rohe, D. 1986, A&A, 159, 276 
Deleuil, M., Meunier, J.C., Moutou, C. , et al. 2009, ApJ, 138, 649 

Dziembowski, W.A., Daszynska-Daszkiewicz L„ Pamyatnykh, A.A. 2007, MNRAS 374, 248 
Fremat, Y., Zorec, J., Hubert, A.-M., Floquet, M. 2005, A&A, 440,305 
Floquet, M„ Hubert, A.-M., Janot-Pacheco, E., et al 1992, A&A, 264, 177 
Foster G. 1995, ApJ, 109, 1889 
Foster G. 1996, ApJ, 111,541 

Ghil, M., Allen, M.R., Dettinge, M.D. 2002, Reviews of Geophysics, 40, 1 
Huang, W., Gies, D R. 2006, ApJ, 648, 451 

Kennely, E.J., Walker, G.A.H., Merryfield, W.J. 1992, ApJ, 400, L71 

Kohoutek, L., Wehmeyer, R. 1999, A&AS, 134, 255 

Neiner, C, Hubert, A. M., Floquet, M, et al. 2002, A&A, 388, 899 

Porri, A., Stalio, R. 1988, A&ASS, 75, 371 

Porter, J.M.,& Rivinius, T. 2003, PASP, 115, 1153 

Rivinius, Th., Baade, D., Stefl, S., et al. 1998, A&A, 333,125 

Rivinius, Th., Baade, D., Stefl, S., Maintz, M. 2001, A&A, 369, 1058 

Smith, MA. 1977, ApJ, 215, 574 

Robertson, T.H., Jordan, T.M. 1989, AJ, 98, 1354 

Steele, I.A., Negueruela, I., Clarck, J.S. 1999, A&AS 137, 147 

Vogt, S.S., &Penrod G.D. 1983, ApJ, 275,661 

Walker, G.A.H., Kuschnig, R., Mathews, J.M., et al. 2005, ApJ, 635, L77 
Zorec, J., Fremat, Y, Cidale, L. 2005, A&A, 441, 235 



